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ABSTRACT 

Recent numerical simulations show that many Lya absorption lines of column 
densities Ahi 10^® cm“^ are produced in transient, mini pancakes. Such pancakes are 
modeled here, approximating the initial perturbation leading to the formation of the 
pancake as a single sinusoidal wave. The density and temperature profiles of the gas 
in the pancake are determined in the case where cooling and heating rate of the gas is 
larger than the Hubble time, which is shown to hold for Zc ~ 3, where Zc is the collapse 
redshift. The Lya absorption line profiles for a line of sight through the pancake are 
then calculated. The absorption lines in general have wings signifying bulk motions in 
the gas. It is shown that the deviation from a single Voigt profile is large for small 
H I column density lines, in which the effect of bulk motions is large. For lines with 
Ahi > 10^^ cm“^, high temperature tend to wash out the signatures of bulk motion. 

The analytical modeling of mini pancakes associated with Lya forest lines— with 
10^^ Ahi 10^® cm“^—gives the corresponding mass scales. In structure formation 
models with CDM and with reheated IGM, the gas in the IGM is inhibited from 
falling into dark matter potential wells with very small velocity dispersions. Structures 
with masses larger than this limit (the Jeans mass) can have infall of gas and form 
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mini pancakes. It is shown here that, for typical values of cosmological parameters, 
absorption lines with N^i ~ 10^^ cm“^ correspond to structures with baryonic mass of 
Mfe ~ Mq with an overdensity of ~ 10 at 2 : ~ 3. The value of A^hi can change by a 
factor ~ 3 in the course of evolution of the pancake in time. It is also shown that there 
is an upper limit to A^hi from a pancake due to the slow recombination rate and the 
importance of collisional ionization at high temperatures. Mini pancakes do not give 
rise to Lya lines with ^ cm“^, for a temperature of the IGM of T = 10^ K, 

J —21 — 1 and ^igm ~ 0.03. 

Subject headings: intergalactic medium - quasars ; absorption lines - cosmology ; 
miscellaneous 


1. INTRODUCTION 

Several ideas for the physical condition of the regions producing QSO absorption lines with 
low H I column density (the Lya forest lines) have been proposed in the past. They have been, 
to name a few, pressure confined gas (Sargent et al.l980; Ostriker and Ikeuchi 1983), gas behind 
intergalactic shock waves (Vishniac and Bust 1987), mini halos (Rees 1986), caustics in velocity 
space (McGill 1990), and other scenarios. Recent numerical simulations have, however, shown that 
the physical condition is a mixture of many such scenarios (Cen et al.l994, Hernquist et al.l995, 
Miralda-Escude et al.l995, Zhang et al.l995). The simulations, performed with a variety of 
structure formation models, show that transient, pancake-like regions of collapsing objects,— and 
not necessarily virialized objects—produce most of the Lya forest lines with Vhi ~ 10^^ cm“^. 
For example, Miralda-Escude et al.(1995) find that the pancake associated with a typical Lya 
forest line has a baryonic mass 10^° Mq at Z 3, with a thickness r\j 50 — 100 h ^ kpc, and 

a transverse size ~ 1 Mpc (where H = 100 h km s“^ Mpc“^ is the Hubble constant at present 
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epoch). The numerical simulations find excellent fits with observations of the H I column density 
distribution of Lya lines and the line profiles. (In the simulations, the Lyman limit and damped 
Lya systems in general come from lines of sight through denser and more spherical objects.) 

That such objects should be prevalent in the universe is not surprising. ‘Zeldovich pancake’ 
is the inevitable result of the fact that a near simultaneous collapse of an object along all three 
directions, as in a spherical collapse model, is very unlikely. Zeldovich (1970) pointed out that a 
one dimensional collapse is more probable to occur. Zeldovich and Sunyaev (1972) then considered 
the fate of gas inside such a nearly two dimensional pancake, although they considered very 
large masses (~ 10^^ M 0 ) which are the typical masses of the first collapsing objects in a HDM 
universe. For a CDM like structure formation model, the first objects to collapse have much 
smaller masses. Therefore, in such a universe the existence of pancakes with smaller masses is not 
unexpected. 

Such pancakes would have masses several times that of the Jeans mass of the intergalactic 
medium (IGM). In a photoionized IGM, with temperatures T ~ 10^ K, the sound velocity is 
~ 15 km s“^. The baryonic gas, therefore, does not fall into, and cool to form galaxies in, dark 
matter potential wells with velocity dispersions smaller than this sound velocity. Recent numerical 
simulations support this scenario (Thoul and Weinberg 1995). The Jeans mass of the IGM, 
therefore, sets the mass scale above which the baryonic gas can fall into the dark matter potential 
well. The Jeans (total) mass is ~ 5 x 10^° h~^ (T/IO^)^'® Mq at 3 for such temperatures in 
the IGM in a Ho = 1 universe. This corresponds to a Jeans length of 0.7 (T/lO^)^'^ h~^ (1 + 

Mpc. This is basically the length Ri c* {1/{Gpcdm)^'^), where the term inside brackets is the 
dynamical time of the halo ( see, e.g., eqn (5) of Reisenegger and Miralda-Escude 1995). The 
(total) masses of the typical first collapsing objects at that redshift would, therefore, be several 
times this mass scale. It is no wonder then that the baryonic masses of the typical mini pancakes 
in the numerical simulations are found to be of the order of 10^*^ Mq. 
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In CDM like theories, smaller structures will also form, but due to the non-zero sound velocity 
of the baryonic gas, the gas will not be clumped in such smaller structures. 

The idea of mini pancakes in the context of Lya lines is also motivated by the recent 
observations of coincident absorption lines in close quasar pairs. The observed fraction (of 
~ 50 — 80%) of coincident lines of two quasar pairs Q0107-025AB and Q1343-I-266AB, with 
proper separations of 360/i“^ kpc and 40/i“^ kpc at redshifts z ~ 1 and z ~ 2 (Dinshaw et al.l994, 
Dinshaw et ah 1995, Bechtold et ah 1994), show that the Lya absorbing systems have transverse 
sizes of the order of a Mpc. With such large transverse sizes, the systems need to have a small 
‘thickness’, i.e. the systems need to be flat, so that the total amount of gas in Lya absorbers do 
not exceed the limit on Of, from Big Bang nucleosynthesis (Rauch and Haehnelt 1995). 

Recently Rees (1995) has discussed the possible mass scales of the pancakes associated with 
Lya forest lines. He has suggested that mass scales between 10® and 10^^ Mq, with various 
collapse factors (in the range of 1 — 3) are enough to cover the range of H I column density of Lya 
clouds below those of Lyman limit systems. The present work is aimed towards quantifying these 
speculations. 

Zeldovich and Sunyaev (1972) discussed the behavior of gas in pancakes approximating the 
initial perturbation that leads to the formation of the pancake as a plane sinusoidal wave. They 
described the density and temperature profiles of the gas as functions of time. The goal of this 
work is to pursue the mini pancake model of the Lya forest lines, along the lines followed by 
Zeldovich and Sunyaev (1972) for large pancakes. Such an analytical picture will give us the 
dependences of the various observables on different parameters, which is sometimes difficult to 
understand from the results of time consuming numerical simulations. The questions asked in this 
paper relate to the physical conditions in these pancakes, the relation between the mass scale of 
pancakes and their H I column density, and also the evolution of the pancake gas in short time 



- 5 - 


scales. 

The plan of the paper is as follows. §2 discusses the physical conditions of the pancake gas 
and determines the density and temperature profile. §3 then calculates the Lya absorption line 
profile and discusses the possible signatures of the pancakes in the observed profiles. The different 
time scales for the pancake gas are then computed in §4, and some of the implications of these 
results are discussed in 55. 


2. GAS IN MINI PANCAKES 
2.1. Preliminaries 

If the mass of the collapsing object is larger than the current Jeans mass, then the baryonic 
gas will also collapse with the dark matter. In the Zeldovich approximation, the trajectory of a 
particle is given by 

r = a(t)(q + 6 (t)p(q)) . (1) 

Here r{t, q) is the actual position of a particle as a function of its Lagrangian (initial) coordinate 
q and time t. The function p(q) describes the initial perturbation, a{t)q is due to cosmological 
expansion, and a{t)b{t) describes the growth of the perturbation in the linear regime. It then 
follows that, in the linear regime, the growth of the density perturbation in three spatial directions 
near a particle with a given q will be given by 

__ 

det((5ij + b{t){dpj/dqi)) 

_pW_ 

[1 - 6 (t)Ai(q)][l - 6 (t)A 2 (q)][l - 6 (t)A 3 (q)] ’ 
b{t){Xi + A 2 + A 3 ) . ( 2 ) 

Here, p{t) is the background density of the universe at t, 5 is the density contrast, and 
—Ai, —A 2 , —A 3 are the eigenvalues of (dpj/dqi). In the generic case, the three eigenvalues will 


P(r,t) 

6{r,t) 
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be different from one another. The probability distribution of the eigenvalues was first derived 
by Doroshkevich. The probability that Ai,A 2 ,A 3 will all simultaneously be either positive or all 
negative is ~ 0.08. The probability that there is contraction in one direction and expansion in 
other two, is 0.42. 

In the generic case, when the eigenvalues are different from one another, the contraction will 
be strongest in the direction for which the corresponding eigenvalue is the largest among the 
three. Eqns (2) show that, the density of the matter around the Lagrangian coordinate q will 
diverge and the matter will be compressed to a sheet in the (actual) Eulerian space, when one of 
the eigenvalues, say, Ai, satisfies b(t)Xi = 1. It can be shown that the density profile in the case of 
one dimensional contraction is given by p oc (see, e.g., eq. (2.7) in Zeldovich and Shandarin 

1989). 


2.2. Gas infall and shock 

Although we use the Zeldovich-Sunyaev formalism for mini pancakes, we must note the 
caveats involved. As mentioned above, mini pancakes are structures with masses just above 
the Jeans scale and, therefore, gas pressure plays an important role, unlike in the case of large 
pancakes discussed by ZS72. However, for analytical simplicity, we neglect gas pressure in this 
work. 

For analytical simplicity, it is assumed that fl = 1. It is also assumed that the gas in the 
IGM is reionized (and reheated) by the time the mini pancakes form. The intensity of the UV 
background radiation (at 912 A) is written as J = J_ 2 i 10“^^ erg s“^ cm“^ sr“^ Hz“^. The 
temperature, T, of the IGM gas roughly obeys T oc (1 + z) after reheating (Miralda-Escude and 
Rees 1994). 

The baryonic gas particles will initially follow the dark matter and form a sheet. A shock wave 
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will soon form which will travel outward. ZS72 (see also Zeldovich and Novikov 1983; Zeldovich 
and Shandarin 1989) estimated the velocity of the shock by assuming that the initial perturbation 
is a plane sinusoidal wave. In reality, of course, the perturbation will be more complicated. 
However, as shown below, the plane wave approximation can be very useful in learning about the 
nature of the gaseous pancakes. 

Pancakes are formed when there is a region in space where the velocity field is coherent. 
Consider a pancake which is associated with such a region with a length scale of perturbation L, 
and therefore, with some mass M (the total mass initially contained in the region of perturbation). 
The perturbation is assumed to be a plane wave with a wavelength equal to A. (Eqn (A2) in 
Appendix A shows the relation between the mass and the wavelength.) The derivation of the 
structure of the pancake using this approximation by ZS72 is briefly described in the Appendix 
A below (after putting the necessary dependences on h and the baryonic fraction of matter H;,; 
note that their results are for a = H = 1 universe). Since the density of the IGM (H/) is more 
relevant for the gas inside the pancakes, than the general fraction of baryons in the universe {0,b), 
the former is used in the equations below. 

Defining a parameter fj, to be the fraction of the total mass M that is inside the region 
bounded by the shock layers, one can write the velocity with which the gas runs into the matter 
in the pancake (in the laboratory frame) (eqn (1) in ZS72; (A7) below), as, 

~ 73.3 (1 + Ws , « 1. (3) 

Here Zc is the redshift of collapse of the perturbation. It can be shown that (for Ho = 1, ZS72; 
(A5) below) 

(.) 3 / 3 ^(^). 2 ^(^)- 2 / 2 ^ ( 4 ) 

where tc corresponds to Zc- For ^ ^ 1, = ^{(t/tc) — 1)^^^. The curve for eqn (4) and the above 



approximation are plotted in Fig. 1 for Zc = 3 and Zc = 4. (It is shown below in §4 that the 
gravitational time scale of mini pancakes limit the possible values of to be ^ 0.5 .) Therefore, fi 
can be used as a parameter equivalent to time t. 


The velocity of the gas inside the shock front is assumed to be zero here, as it is expected to 
be small (see Bond et al.l984). The velocity of the shock front is then given by, Vgh = (1/3)17 
(Jones et al.l981. Bond et al.l984). Note that U ss Vgh in ZS72. 


For small values of jj,, the (particle) gas density in front of the shock is (eqn (3) in ZS72; (A6) 
below), 

{l+Zcf 


m 


3n, 


O 9 9 


// < 1 , 


= 2.3x10-U-2(1 + z,)3/x-2(^) 


(5) 


where it is assumed that the IGM gas is ionized {p = O.Qnmp) and Uq is the particle density 
at present epoch. The postshock gas density is Ugh = 4ni, for a strong shock. The postshock 
temperature is, again for small p, 


Tsh = 3 ^ = 1.3 X 10® (1 + Zc) ( 


M ^2/3 

p K 


IQii Mq 


( 6 ) 


The pressure behind the shock front {psh = (4/3) pi t/^) is a independent of p (i.e., of how much 
of matter is inside the shocked layer for a given pancake ) for small p (ZS72). 


The distance of the shock front from the central plane of the pancake can be estimated as. 


Tsh 


= J Vshdt = J Vsh ^ dp 

« 0.27/1-2/3 + 


p 1 . 


(7) 


As mentioned above and explained in the Appendix, the above expressions are obtained in 
the absence of gas pressure. As gas pressure is more important for mini pancakes than the large 
pancakes studied by ZS72, these expressions are only very approximate. 
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2.3. Density and temperatnre profile 


What is the maximum gas density in the pancake? The maximum density would be infinite 
for infinite cooling rate and for a zero pressure gas outside, neither of which is true in the present 
case. It can be shown that the cooling time scale of the gas in the presence of the UV background 
is larger than the Hubble time at 2 : ~ 3 for the relevant gas densities of this problem (see Appendix 
B). 


However, the heating time scale becomes comparable to or shorter than the Hubble time 
for Zc ~ 3 (Appendix B). In the following, we restrict our analysis to Zc ~ 3 where both cooling 
and heating time scales are larger than the Hubble time and assume that the gas contracts 
adiabatically. 


Furthermore, since the evolution of the IGM temperature T depends on various parameters 
and models of ionization (Miralda-Escude and Rees 1994), we will, for simplicity, assume a 
constant IGM temperature of T ~ 10^ K in the following. Introducing an evolution of T will 
introduce more free parameters, which we will avoid. 


If the ambient IGM density is n, with pressure p (and temperature T), then the maximum 
gas density and temperature in the pancake is given by 


Ur 


= T 


fPsh\^/^ _ 

1.2xl0-^hi2/5j 

T +ZcN^18/5 

r ^ ) 

\ p ) 

. 4 ) ' 


f IT'max 

= 4.7 X 10^ 

f ^ 4- Zc\ 2/5 j 

r M N 

V n ) 

V 4 / ' 



2/5 / Qj ^ 

votT^ 

n4/15 / f \ 

) Vio-^K? 


' T n-3/5 

• lO^K? 

3/5 


( 8 ) 


As mentioned in §2.1, in the Zeldovich approximation, the density of the gas has a profile 
(in the direction perpendicular to the pancake) that is given by p oc (e.g., eqn. (2.7) of 

Zeldovich and Shandarin 1989). The density profile of the gas can, therefore, be divided into four 
characteristic regions: (a) in front of the shock, the density follows a profile, as explained 

above; (b) the density jumps by a factor of 4 at the shock front; (c) behind the shock, the gas 
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again follows a profile, in the Zeldovich approximation, as explained above; (d) and then it 

reaches a maximum when the density is Umax of eqn (8). For analytical simplicity, we assume this 
central region, say, r < Vq, with n ~ rimax to be a region of uniform density. 

In the single wave approximation used here, we have neglected the motion of the gas in the 
directions along the pancake. In reality, gas will also move along the pancake. Modelling this 
effect is, however, beyond the scope of the single wave approximation. We neglect it here, and 
assume that the density in the central region is a constant in time, for the short time evolution we 
are interested in here (// <C 1). Gas will also tend to expand in the perpendicular direction, but 
for small n (which is considered in the present work), pressure behind the shockp^/i, stays constant 
with fj, (that is, in time) and the ram pressure will tend to confine the gas. In other words, we 
assume that. 


n ~ 

'^max 


r <ro 

oc 

, r , 
f^max f , 
To 

)- 2/3 

Vo <r < rgh 


'^sh 


r ~ Tgh 

oc 

Ugh , r 

4 

- 2/3 

r> rgh- 


(9) 

If rsh < fo, then n ~ Umax, r < ^sh, and n oc ^ ^ ^ Matching the densities at 

the various boundaries, as explained above, one obtains for the extent of this central region, as 

- 4/15 


+ Zc\ 

- 19/10 f M 

V 4 ) 

UOIIMq/ 


kpc . 


( 10 ) 


This length scale will be important in calculating the column densities in §2.4. 


Similarly, the temperature has a maximum value in the central region, as given in eqn (9). At 
the current position of the shock, Vsh, the temperature will be given by eqn (6). For Tq < r < Vsh, 
the gas would have been shocked and would be raised to corresponding Tg^. As Tg^ oc and 
rgh oc it is easy to see that T oc for Vq < r < Vgh- Outside rgh, the temperature the same as 



-11 - 


that of the IGM (T). We remark here that the temperature profile will in general change in time 
due to cooling. As explained above, the cooling rate is small for z 3 and the above temperature 
and density profiles hold only in that regime. 

The profiles of (5 = n/h, 5hi = nni/nHi and T/T for a typical pancake are shown in Fig. 

2, for jj, = 0.3,0.4. We have used M = 5 x 10^^ Mq, h = 0.75, Zc = 3 and 14/ = 0.03. For each 
value of iJ,, the values of and rgh are first evaluated (eqns (7) and (10)), and then the density at 
each point in space using eqns (5) and (8). The density of neutral hydrogen is evaluated at each 
point of the pancake using ionization balance (see §2.4). The increase in the overdensity at later 
times is due to the decrease in the ambient density in the universe in time. It is seen that the gas 
in the pancakes have a typical maximum overdensity of order ~ 10-50 (see also Miralda-Escude 
et al.l995). One should note that in reality gas will expand along the pancake, thereby decreasing 
the central density. Therefore, the overdensity obtained above should be taken as an upper limit. 

The density and temperature profiles in Fig. 2 are much different from those in 
large, massive pancakes, as studied by ZS72, in that the maximum density in their case is 
much larger, and therefore, Tq is smaller (it was small enough to be unimportant for those 
pancakes). Also, T^ax in ZS72 was much lower because a cold IGM gas was assumed. 

For small fi, the temperature profile does not have the usual ‘double horn’ profile that is 
expected in large pancakes (Jones, Palmer and Wyse 1981). Such a profile is possible if, 
fi ^ 0.3((1 + (M/10^^ (T/lO^ As discussed below, large values of 

[I in general correspond to large A^hi- Therefore, such temperature profiles are not expected for 
small H I column density Lya lines. The spectra and the corresponding density and temperature 
profiles obtained in the numerical simulations of Miralda-Escude et al.(1995) support this. 

The extent, or the thickness, of the pancake can be estimated by finding the distance r' where 
the density contrast 5 = {n — n)/n (where n is the background IGM density) is unity. Eor the 
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above profile, this happens at r' ~ rgh ~ 6.2 x 10 ^ ^ For M = 3 x 10^^ Mq 

(corresponding to M;, 10^° M©, for VLi 0.03), h = 0.5 and Zc = 3, the ‘thickness’ of the pancake 

is ~ 20 kpc. For comparison, the typical thickness for pancakes as reported by Miralda-Escude 
et al.(1995) is 50 — 100 h~^ kpc. It shonld be noted that this is also consistent with the upper 
limit on the thickness (~ 50 kpc) put from the combined knowledge of the large transverse sizes of 
Lya absorption systems and the value of from Big Bang nucleosynthesis (Haehnelt and Rauch 
1995). 


2.4. Neutral and H I column density 


For a line of sight which is perpendicular to the pancake, the total column density of the gas 


is. 


rr=rsh 


~ 2 / n{r)dr =2 Umax To {1 + 8.8 ^ fi} 


/r=0 


( ^ \ ( h \4/3 .l + ZcN 

; / M s 

.1/3 / c 

V0.3/ V0.5/ V 4 / 

VlOiiM©^ 

' Vo.oJ 


cm 


o 

-2 


( 11 ) 


This is comparable to the values of column density plotted in Fig. 9a of Miralda-Escude 
et al.(1995). Eor a gas in ionization equilibrium with photoionization by the UV background 
radiation, and with collisional ionization, one can also estimate the H I column density. 
Photoionization and recombination rates are taken from Black (1981), and Cen (1992). Fig. 3 (a) 
and (b) plot N^i and Nh vs. the mass (total and baryonic) that corresponds to the pancake, for a 
line of sight perpendicular to the pancake. Here Mi, = fj,Q,j M. Fig 3 (c) plots Nh against N^i. 
Fig. 3 (d) shows the evolution of in time, or equivalently, with increasing values of fi. 


The plots show that, in general, larger (more massive) pancakes correspond to higher IVhi 
( approximately, N^i oc (see below in §4 for a discussion on this). There is however a scatter, 

due to possible differences in the ambient density (H/), collapsing redshift (zc), and due to the 
evolution of pancakes in time. Since Nh oc Qj, for a photoionized gas, N^^ oc (see below). 
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The column density increases at higher collapse redshifts, but this can be an artifact of the 
assumption that a constant background temperature (T). In reality, T may increase with redshift 
(in a photoionized universe T oc (1 + z) much after the reheating epoch). This will decrease the 
density in the central region and will decrease the column density. Such an evolution in the IGM 
temperature introduces another parameter for the pancakes, and will not be discussed here to 
avoid confusion. The evolution of the individual pancakes first increases A^hi till about ^ ~ 3 
(for M ^ 10^^ Mq), after which hotter temperatures deplete more neutral atoms and halts the 
increase in 

This behavior is also reflected in Fig. 3(c) which plots Nh against A^hi- Large values of A^hi 
are either associated with large M or large fx, both of which are associated with high temperature, 
which tends to decrease the neutral fraction. This fact, therefore, puts an upper limit on the value 
of A^hi for Lya lines from pancakes. In particular, for J_ 2 i = I and 0,j ~ 0.03, mini pancakes can 
only give rise to Lya lines with A^hi ~ 10^^ cm“^ for Zc = 3. For Zc = 4, this limiting value is 
found to be A^hi ~ 10 ^^'^ cm“^. 

3. ABSORPTION LINES 

Armed with the knowledge of the density and temperature profiles, one can calculate the 
Lya absorption line profile for a line of sight through the pancake. The additional information 
necessary is that of the velocity field. The bulk velocity of the gas within the shocked layers 
is small and is neglected here. Outside the shocked layer, the gas particles stream towards 
the pancake with velocity that is determined by the initial velocity perturbation. In Zeldovich 
approximation, the velocity of a particle remains constant. This allows one to find the velocity 
profile of the gas at later times. 


Consider a particle initially at Lagrangian (comoving) coordinate ^ with velocity v. This 
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particle will eventually meet the shock front with the same velocity v at some point r in the 
Eulerian (physical) coordinate. This velocity v and the point r are given by eqns (3) and (7). Also, 
as described in Appendix, every point ^ is associated with a value of Since the velocity remains 
constant, this mapping is also true for other points in physical space. For any other point r', one 
needs to find the corresponding or from eqn (7) and then find the velocity from eqn (3). 

Let r = 0 be the origin (in the dimension perpendicular to the pancake) at some redshift z. 
The total velocity of a particle at r is therefore, V (r) = Hr + v{r) where H is the Hubble constant 
at z and v{r) is the peculiar velocity. The profile of the absorption line is determined by the 
overdensity of H I in velocity space, convolved with a Gaussian of width b{r) = {2kT{r)/mp}^'^. 
The smoothed H I overdensity at some velocity V in velocity space is given by (Miralda-Escude 
and Rees 1992), 

<^hi(R') = J dr5ai{r)H exp(- ^ ^ bir)) . (12) 

If the Gunn-Peterson optical depth to Lya absorption by the IGM is tgp, then the total optical 
depth at velocity V is given by t{V') = tgp{^ + (5hi(R 0)- Here, tgp is the Gunn-Peterson opacity, 
and in a photoionized universe it is given by tgp = 0.32 (1 + z)*'^ ( using eqn (39) of 

Miralda-Escude and Ostriker (1990), with clumping factor / = 1 and the spectral index a = 1). 

Examples of line profiles are plotted in Fig. 4, for various values of fx for a particular pancake, 
and for different orientations. Also shown are the fits to the profiles using Voigt function, which 
depend on and a b parameter. The value of Vhj changes when the angle of inclination (0)of 
the pancake to the line of sight is changed, and is roughly proportional to (cos(0)). (The total H I 
column density is exactly proportional to (cos(0)), but this is only approximately true for the H I 
column density of the fitted central line, as some of the neutral atoms also show up in the extra 
lines.) This is mainly because of our assumed flat density profile for r ^ ro¬ 


ll is seen that the absorption lines from pancake can be fit by Voigt profiles, although with 
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some deviations. In general, the absorption profile consists of a central dip, which follows a Voigt 
profile, and some extra dips which are shifted from the central one in frequency space. These wings 
make the value of Vhi from fitting the central portion of the profile smaller than the corresponding 
value from Fig. 3 (d). Absorptions by some neutral atoms show up in the profiles away from the 
central portion - causing the absorption profile to deviate from the Voigt profile. 

Even though pancakes are associated with bulk motions, the thermal gas makes the profile 
close to Voigt profile for lines with large iVni. Lines with small iVjji, however, are in general 
associated with small values of when only a small fraction of the total mass has been heated 
through adiabatic compression. One, therefore, expects to see the deviations from a simple Voigt 
profile in the case of low mass pancakes, or when ^ is small in a pancake, when bulk motion 
dominates over the thermal contribution to b. Such an example is shown in Fig. 4(d). It is similar 
to some of the low column density lines in the spectra obtained in the numerical simulations of 
Miralda-Escude et al.(I995) (see their discussion in §3.3.3). 

The deviation of the absorption profile from a Voigt profile can be quantified by the ratio of 
the equivalent width of the central dip with the Voigt profile (IFc) to the total equivalent width 
in the absorption profile (Wtot)- As shown in Eig. 5a, this ratio is very small for lines with small 
total equivalent widths. These lines are produced by less massive pancakes and for pancakes with 
small values of ^ (and therefore with smaller values of Vhi). The curves show that a fraction 
of O.I of the total equivalent width of the absorption line shows up in the wings. The wings 
could, however, be fitted separately as extra lines with Voigt profiles, but the interpretation of 
these (small N^i) lines as separate ‘objects’ or ‘clouds’ would not be appropriate, as these extra 
lines simply reflect the bulk motion in the pancake (see also Rauch 1995). For M ^ 5 x 10^^ Mq, 
the ratio Wc/Wtot ~ 0.5 for Wtot ~ O.I A, corresponding to iVni 2 x 10^^ cm“^. Several of the 
observed small Vhi (IVhi ~ 10 ^^ cm“^) lines, therefore, reflect the bulk motions relevant for the 
formation of mini pancakes (see also McGill 1990, where it was pointed out that velocity cusps in 
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redshift space can masquerade as Lya lines). 

For small deviations from the Voigt profile, the b value of the central line is determined by 
thermal motion of the atoms. Fig 5b plots the ratio of the b from Voigt profile fitting of the 
central line, to the b that corresponds to Tmax of eqn (6), as a function of Vhi of the central line 
(determined from profile fitting). It clearly shows that for lines with ^ 10^^ cm“^, the b value 
is determined by the thermal motion of the atoms. For lines with ^ 10^^ cm“^, however, the 
temperature in the shock front dictates the value of b (and therefore deviates from the value of 
b corresponding to T^ax of the central region). A statistical study for the distribution of b is, 
however, beyond the scope of the present work. 

4. DISCUSSIONS 

In CDM type scenarios the mass scale of the first gaseous pancakes is given by the Jeans 
mass. The cutoff in the power spectrum, however, is not as sharp as the lower cutoff in HDM 
scenarios. Also, the power spectrum in the vicinity of the Jeans cutoff is flatter for CDM type 
structure formation models than the HDM power spectrum. This means that the masses of the 
first gaseous pancakes may not have a single value and is likely to have a distribution. 

It is not clear whether or not the column density distribution of Lya lines arises from this 
distribution in mass. If that is the case, then the above relation can be used, in principle, to 
calculate the distribution of Ahi (number of lines with H I column density Vhi per unit H I column 
density per unit redshift, or {dN/dN^n dz)), if an expression for the distribution function of the 
pancakes can be found in terms of their masses (or, equivalently, the length scale L of the region in 
which the initial velocity field is coherent). However, the cause of the column density distribution 
of Lya lines is likely to be more complicated than this. Firstly, if the pancakes survive for about 
a Hubble time, the evolution of the pancakes need to be taken into account. In any case, such a 
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statistical study is beyond the scope of this work. 

Several observations involving double lines of sight suggest that the sizes of regions in the 
direction perpendicular to the line of sight can be several hundreds of kpc. Miralda-Escude 
et al.(1995) find the typical transverse size of pancakes to be ~ 1 Mpc. As mentioned earlier, the 
Jeans scale at z ~ 3 is Aj ~ 0.13 (T/lO^ K)^-^ (/i/0.75)“^ ((1 + z)/4:)~^-^ Mpc, for a photoionized 
universe with T oc (1 + z). The typical first pancakes with infalling baryonic gas will therefore have 
sizes several times this (the diameter will be twice this value). It is, therefore, conceivable that the 
transverse size of pancakes would be ~ 1 Mpc. More over, the gas flows in the transverse direction 
will make the size even larger. As mentioned earlier, when one of the eigenvalues (Ai, A 2 , A 3 ) is 
negative (which describes contraction), the probability that the other two eigenvalues are positive 
(that is gas expands in the transverse direction), is ~ 0.42 (ZS72). As suggested by many previous 
authors, this expansion can make the pancakes larger in the transverse direction, and increase the 
probability of detection. 

Finally, a word about the approximation used is in order. The results obtained here should 
be considered as only very approximate, and have large error bars, because of the approximation 
of single plane waves. However, it is interesting that some information about the behaviour of the 
gas in pancakes, and the effects on the absorption profile, can still be obtained from such simple 
approximations. 


5. CONCLUSIONS 

In this paper, mini pancakes are analytically modeled with single plane wave perturbations, 
and the consequences for the Lya absorption lines from these mini pancakes are discussed. 
Although a very crude approximation to the real pancakes, this model is shown to elucidate a 
number of important points. The most important results of the paper are contained in Figs. 3 
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and 5. They show the mass scales of pancakes that corresponds to the HI column densities of Lya 
forest lines with A^hi ~ cm“^, the evolution of A^hi in time for different pancakes, and the 

importance of bulk motions in low equivalent width absorption lines. It is shown that pancakes 
with total masses in the range of M ~ 10^^-1012 . Mq (and a baryonic mass of Mi, ~ 

Mq ) can give rise to Lya lines with A^hi ~ cm“^. Here, the mass corresponds to the size 

of the perturbation. These structures have an overdensity of the order of ~ 10. Also, there is a 
maximum value of Ahi for a given value of H/. For J_ 2 i = 1, a temperature of the IGM of T = 10^ 
K, and H/ ~ 0.03, mini pancakes only give rise to Lya lines with A^hi ~ 10^^'^ cm“^. These results 
hold only for Zc ~ 3 for which the gas in the mini pancake can be assumed to be compressed 
adiabatically. 

Calculating the expected absorption profiles from the pancakes, we also confirm that the 
information of the bulk motion associated with pancakes is carried by low HI column density lines 
with Ahi ~ 10^^ cm“^. The lines in general consist of a central dip with wings that signify bulk 
motions. The wings are most significant for low equivalent width lines. The lines could be fitted 
with more than one Voigt profiles, but would not necessarily mean that the lines arise in separate, 
discrete objects. 

I am grateful to my colleagues in lUCAA and Dr. K. Subramanian for valuable discussions, 
encouragement and criticisms. I also thank the anonymous referee whose comments have improved 
the manuscript. I am indebted to Prof. P. J. E. Peebles for an inspiring discussion on Lya lines 
and pancakes. 
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Figure captions Fig. 1 - The parameter fi is plotted against redshift z for two different collapse 

redshifts Zc = 3,4. Solid curve is the result of eqn (4) and the dashed curve corresponds to the 
approximation used in the work. 

Fig. 2 - (a) Curves show the overdensity profile of pancakes with M = 5 x 10^^ Mq at // = 0.3 
(solid curve) and at /x = 0.4 (dotted curve). Most of the increase in the overdensity in time comes 
from the decrease in the ambient density. 

(b) Similar curves for H I overdensity profile. 

(c) Curves for the temperature profile for the same parameters. The temperature of the IGM gas, 
f = 10^ K. 

Fig. 3 - (a) The H I column densities A^hi of pancakes, for regions where the H I overdensity 
is greater than unity, are plotted against the total mass M of the pancakes for lines of sight 
perpendicular to the pancake. Dotted, short dashed and long dashed curves refer to pancakes with 
/X = 0.2,0.3, 0.4 respectively, for D/ = 0.03. The collapse redshift of the pancakes Zc = 3. 

(b) Same as in Fig. 3(a), except that A^hi is plotted against the baryonic masses Mb = fiQjM. 
The upper set of curves correspond to collapse redshift Zc = 4: and the lower set, to Zc = 3. 

(c) Total column density Nh is plotted against A^hi for fi = 0.2, 0.3,0.4, Omegai = 0.03, Zc = 3. 

(d) The H I column densities A^hi of different masses are plotted as functions of g (or, equivalently, 
time). From top to bottom, the masses involved are M = (50,10,5,1) x 10^^ Mq. The dotted 
curve shows the contour left of which which the total equivalent width Wtot of the corresponding 
Lya absorption profile is less than 0.1 A. 
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Fig. 4 - (a) The Lya absorption profile for a pancake with M = 5 x 10^^ Mq, /i = 0.3 for a line of 
sight perpendicular to the pancake, for Zc = 3., h = 0.75, 12/ = 0.03. The Voigt profile fitted to 
the central dip has = 2.4 x 10^^ cm“^ and 6 = 35 km/s. 

(b) The profile for the same pancake but with ^ = 0.4. The Voigt profile has Vhi = 10^^ and 
6 = 35 km/s. 

(c) The profile for the same pancake with ^ = 0.3 but for a line of sight making an angle 45° to 
the pancake. The Voigt profile fitted has Vhi = 3.4 x 10^^ cm“^ and 6 ~ 35 km/s. 

(d) The profile for a pancake with M = 10^^ cm“^ at /U = 0.25 (the cosmological parameters 
remain the same). The Voigt profile fitted to the central dip has N^i = 3.6 x 10^^ cm“^ and 6 = 30 
km/s. 

Fig. 5 - (a) The ratio of the equivalent width of the fit to the central dip Wc to the total equivalent 
width of the absorption profile Wtot is plotted against Wtot for different masses. Curves with long 
dash, short dash and solid line correspond to M = 10^^, 5 x 10^^, 10^^ Mq with h = 0.75, Zc = 3 
and 12/ = 0.03. 

(b) The ratio of the b value of the Voigt fit to the central line to the b corresponding to T^ax of the 
pancake is plotted against the values of corresponding Vhi of the central line. The cosmological 
parameters are the same as in Fig. 5(a). 

Fig. B - The heating and cooling rate of gas with density n// = 10“^ cm“^ in the presence of a 
photoionizing background of intensity J_ 2 i = 1 are plotted as functions of the temperature. The 
upper solid curve is the cooling rate for J_ 2 i = 0. The solid line at bottom is the cooling rate for 
J_ 2 i = 1 and the dashed line is the heating rate for the same ionizing background. 
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A. APPENDIX A 


The derivations of eqns (3), (4), and (5) are sketched briefly below, following ZS72 (see 
also Zeldovich and Novikov 1983). The perturbation leading to the formation of the pancake is 
assumed to be a sinusoidal plane wave in an = 1 universe, with baryon fraction fife. The motion 
of particles in the direction of contraction (before being compressed by shock waves) is given by 

r = _ 5^4/3 gin _ (Al) 


Here ^ is the Lagrangian coordinate. The mass M that is contained in the region (of comoving 
length scale A) of perturbation is given by. 




A(1 + z) \ 3 

-,p = 


1 




;A = 


27rt2/3 

k 


(A 2 ) 


Here, po = 1.89 x 10^® h'^ fl/ is the present day mass density of the IGM. For matter not yet 
compressed by the wave, the density is given by 


1 


1 


^ GvrGt^ 1 — 61^/3 A: cos^ ^ 

2/3 

Therefore, tc, the epoch of the formation of the pancake at r = 0 is determined by kbtc = 1- 
The parameter p, is then defined as p = This means that the interval 0 <^< 7 r/A:in the 
Lagrangian coordinate is equivalent to 0 < /x < 1. 


If the pressure of the IGM gas is assumed to be zero, then a plane at r = 0 reaches infinite 
density. One can think of a wave front in the Lagrangian coordinate corresponding to the particles 
reaching this plane. This wave front moves outward, and the Lagrangian coordinate pi of the 
wave front is given by, 

r = 0 = TTpi — {hkty^) ^ sinyr^i . (A4) 

Note that, in reality the shock will form at a distance from the r = 0 plane, and the above 
expression, which assumes an origin at r = 0 is reasonable only for small gas pressure. However, 
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as explained in the text, we neglect gas pressure in the present work. This shows that, 

t\2/3 / TTfMi 


(Ly^ - ( \ 

\tc^ VsinTT^i/ 


(A5) 


Using (A3), the density at the front is then given by, 
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(A6) 


The velocity of matter running into the wavefront is obtained as, (for small and for t ~ tc, 
i.e., ^ <C 1) as. 
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(A7) 


Here {1 + Zc) corresponds to tc- 


For matter with non-zero pressure, the above equations will remain approximately unchanged 
for large fj,, when the wave has moved a large distance. The central regions will, however, not 
attain infinite density in this case (see §2.3). 
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B. APPENDIX B 

We present the cooling curve for a gas in the presence of ionizing UV photons here. The 
cooling processes used in the calculation include collisional ionization, recombination, dielectronic 
recombination and collisional excitation, in a gas containing hydrogen and helium (X = 0.76), as 
listed in Black (1981). We used the photoionization cross-section for J_ 2 i = 1 as in Black (1981) 
(for a spectral index a = 1). Figure B shows the cooling rate (solid lines) for a gas with nn = 10“^ 
cm—3 with J_ 2 i = 0 and J_ 2 i = 1, and shows that line cooling is suppressed by the ionizing 
photons substantially. For higher densities, the suppression of line cooling is less. Therefore, for 
a maximum density of ~ 10“^ cm“^ and temperature 10® K, as in eqn (8), the cooling rate of 
the gas is ~ ~ 10“^^ erg cm“® s“^. This implies a cooling time scale of 10^^ (T/IO® k) s, 

which is larger than the Hubble time at z ~ 3. 

Fig. B also plots the heating rate (dashed line) for the same density and for J_ 2 i = 1. It is 
seen that at T ~ 10^'^ and n = 10“^ cm“® (eqn (8), the heating time scale is ~ 3 x 10^® s, 

comparable to the Hubble time at z ~ 3 for a flat universe and h = 0.75. 

Therefore, the assumption of adiabaticity for the gas in mini pancakes is reasonable only for 
Zc ~ 3. 
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